Can Polymer Coils be modeled as "Soft Colloids"? 
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We map dilute or semi-dilute solutions of non-intersecting polymer chains onto a fluid of "soft" 
particles interacting via a concentration dependent effective pair potential, by inverting the pair 
distribution function of the centers of mass of the initial polymer chains. A similar inversion is used 
to derive an effective wall-polymer potential; these potentials are combined to successfully reproduce 
the calculated exact depletion interaction induced by non-intersecting polymers between two walls. 
The mapping opens up the possibility of large-scale simulations of polymer solutions in complex 
geometries. 
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A statistical description of polymer solutions in com- 
plex geometries, such as the colloid-polymer mixtures 
which have recently received much experimental atten- 
tion generally relies on a nanometer scale segment 
representation of the polymer coils, a computationally 
very demanding task except in the special case of ideal 
(non-intersecting) polymers obeying Gaussian statistics 
j|. This obviously follows from the fact that, although 
the colloidal particles may reasonably be modeled by 
hard impenetrable spheres or other complex shapes lack- 
ing internal structure, each polymer coil involves L seg- 
ments which must satisfy a non- intersection constraint. 
It thus appears natural to attempt a mesoscale coarse- 
graining, whereby polymer coils interact via effective pair 
potentials acting between their centers of mass (CM). 
Since polymers can interpenetrate, the effective potential 
(3v(r) is expected to be soft, with a range of the order 
of the radius of gyration R g of individual coils. Such 
a coarse-grained description has been a long-time goal 
in the statistical mechanics of polymer solutions, dat- 
ing back to the first attempts by Flory and Krigbaum |^] 
who employed mean-field theory to find an interaction for 
which the strength at overlap scales as: /3v(r=0) ~ L 2 . 
Later, scaling arguments jqL field-theoretical renormal- 
ization group calculations [Q , and simulations S] con- 
firmed that the range of the interaction between two iso- 
lated polymer coils is of order R g , but found that in the 
scaling limit the strength (3v(r=0) is independent of L 
and of order fc^T. 

In this letter, we show that a meaningful "soft colloid" 
picture of polymer coils may be built on a coherent "first 
principles" statistical mechanical foundation. We derive 
both the effective wall-polymer CM interaction (3(f)(z), 
and the "best" local effective pair-potential (3v(r) be- 
tween polymer CM's for finite polymer concentrations. 
These potentials are then applied to simulate bulk poly- 
mer solutions, as well inhomogeneous polymers near a 
hard wall and polymers confined between two parallel 
walls to extract the effective depletion potential between 
plates. The "soft colloid" approach turns out to be suc- 



cessful not only in the dilute regime but also, perhaps 
more surprisingly, well into the semi-dilute regime. A 
related "soft particle" picture has been applied to poly- 
mer melts and blends but the corresponding phe- 
nomenological implementation differs substantially from 
the present "first principles" approach. 

We consider a popular model for polymers in a good 
solvent [ fio|| , namely N excluded volume polymer chains 
of L segments undergoing non-intersecting self avoid- 
ing walks (SAW) on a simple cubic lattice of M sites, 
with periodic boundary conditions. The packing frac- 
tion is equal to the fraction of lattice sites occupied by 
polymer segments, c=N x L/M, while the concentration 
of polymer chains is p=c/L=N/M. For a single SAW 
chain, the radius of gyration R g ~ L u , where v ~ 0.6 
is the Flory exponent p0[ . The overlap concentration 
p* , signaling the onset of the semi-dilute regime, is such 
that Anp*Rp3 ~ 1, and hence p* ~ L~ 3 " . We have 
carried out MC simulations for chains of length L=100 
and L=500, and covered a range of concentrations up 
to pj p* (~ 5. The pair distribution function g(r) of the 
centers of mass was computed for several concentrations; 
g(r=0) is always non-zero, thus confirming the "softness" 
of the effective pair potential (3v(r). The latter was then 
derived from g(r) by an inversion procedure based on the 
hyperneted-chain approximation (HNC) closure relation 
& 



g(r) = exp{-/?u(r) + g(r) - c(r) - 1}, 



(1) 



where 0=1/ ksT, while c(r) is the direct pair correlation 
function, related to g(r) by the Ornstein-Zernike (OZ) 
relation B]. To any given g(r) and density there corre- 
sponds a unique effective pair potential /3v(r), capable of 
reproducing the input g(r), irrespective of the underlying 
many-body interactions in the system Jl2| ; in a variational 
sense this f3v(r) provides the "best" pair representation of 
the true interactions jl3| , and leads back to the true ther- 
modynamics via the compressibility relation pd[ | . While 
the simple HNC inversion procedure would be inadequate 
for dense fluids of hard core particles, where more so- 
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phisticated closures or iterative procedures are required 
p3| , we are able to demonstrate the consistency of the 
HNC inversion in the present case [Q. If the resulting 
effective (3v(r), examples of which are shown in Fig. [l], 
are used directly in MC simulations, the calculated "ex- 
act" g(r) for this effective representation coincides within 
statistical errors with the g(r) derived from the simula- 
tion of the full initial polymer segment model. In fact, 
the HNC closure turns out to be quasi-exact when ap- 
plied to the simple Gaussian model ]l5j whereby parti- 
cles interact via the potential f3v{r)—eexp[~ a (r/R g ) 2 ], 
which yields a reasonable fit to the effective pair poten- 
tials shown in Fig. [l]. Even the much cruder random 
phase approximation closure, c(r)=— (3v(r), yields semi- 
quantitatively accurate results in the regime of interest 
fn3| , fl7f . Careful inspection of Fig. [j] reveals that the ef- 
fective pair potential is not very sensitive to the polymer 
concentration. The value at r—0 first increases slightly 
with p, before decreasing again at the highest concentra- 
tion. More strikingly, and perhaps not surprisingly, the 
range of (3v(r) increases with p. The effective potential 
becomes slightly negative (0(lO _3 fcBT)) for r/R g > 3 at 
the higher concentrations. 
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FIG. 1. The effective polymer CM pair potential 
j3v(r/R g ) derived from an HNC inversion of g(r/R g ) for 
different densities. The x-axis denotes r/R g , where R g 
is the radius of gyration of an isolated SAW polymer. 
Inset: The polymer CM pair distribution function g(r) 
calculated for L=500 SAW polymers and used to gener- 
ate /3v(r). 

The properties of soft-core fluids are significantly dif- 
ferent from their hard-core counterparts. For example, 
for potentials of the type shown in Fig. (Q), the pressure is 
very well described by (3P=p+l/2(3V(0)p 2 over the entire 
density range |T^,[l7j]. Here V(0) is the Fourier transform 
of the potential, at k=0. Our observation that potentials 
become slightly longer ranged at higher densities implies 
that the pressure scales with an exponent slightly higher 
than 2, so that the equation of state (e.o.s.) is at consis- 



tent with the well-known p 225 law p0| . At first sight it 
may seem surprising that a two-body potential could re- 
produce the full e.o.s. without explicit many body terms. 
However, the effective potential we use is constructed to 
reproduce the true thermodynamics through the com- 
pressibility relation (ignoring small volume terms); the 
relative insensitivity of (3v(r) to concentration implies 
that many-body interactions are not very important |l6[ . 

This insensitivity to concentration makes it possible to 
apply the effective potential appropriate for a given mean 
concentration to inhomogeneous cases, where the local 
polymer concentration deviates from the mean. Such a 
situation occurs when a polymer solution is confined by a 
hard wall. Using the same explicit SAW polymer model 
in MC simulations, we have computed the "exact" pro- 
files h{z)=p{z) I p — I, where z denotes the perpendicular 
distance of the polymer CM from the wall. 




FIG. 2. The wall-polymer potential f3<f}{z/R g ) de- 
rived from an HNC inversion of h(z). Inset: The 
wall-polymer density profile h(z)=p(z)/ p — 1 for dif- 
ferent densities. The corresponding adsorptions V are 
0,0.096,0.132, 0.178,and 0.248 in units of R g 2 respec- 
tively. 

Examples of h(z) for several bulk concentrations are 
shown in the inset of Fig. |2[ The corresponding adsorp- 
tions r are defined by: 

Jo 

where VL ex /A is the excess grand potential per unit area, 
p the bulk concentration of the polymers and p their 
chemical potential. From a knowledge of the concentra- 
tion profile p(z), and the bulk direct correlation func- 
tion between polymers CM's c(r), one may extract an 
effective wall-polymer potential f3<j){z) by combining the 
wall-polymer OZ relations jll) with the HNC closure, re- 
sulting in: 

= ^ MF (z) +p [ dr'h(z')c{\v - r'|). (3) 
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The first term is the usual potential of mean force 
@4> MF (z)=— In [p(z)/p], to which f3<j){z) would reduce in 
the p — » limit, while the second term arises from corre- 
lations between the polymer coils next to the wall. Us- 
ing the c(r) extracted from the earlier bulk simulations 
of g(r), together with Eq. (||), we are able to extract 
(3(j>(z) from the density profiles. Results for various bulk 
concentrations are plotted in the Fig. (0). 

The range of the effective wall-polymer repulsion in- 
creases with increasing concentration, while the density 
profiles actually move in closer to the wall. It is impor- 
tant to stress that the correlation term considerably en- 
hances the repulsion compared to the potential of mean 
force. We have tested the consistency of the inversion 
procedure (which, to the best of our knowledge, has not 
been attempted before for any wall/fluid interface) by 
using f3<f>(z), and the pair potential (3v(r) for the appro- 
priate bulk concentration, in MC simulations based on 
these effective interactions (such simulations are at least 
an order of magnitude faster than simulations of the ini- 
tial segment model). The resulting concentration profile 
of the effective "soft colloids" agrees to within statistical 
accuracy with the initial p(z) obtained from the detailed 
segment simulations, and the corresponding adsorption T 
differs by less than 1% from the exact value, thus demon- 
strating the adequacy of the "soft colloid" representation 
of the interacting polymer coils. 

An even more severe test of this representation is pro- 
vided by a calculation of the depletion interaction be- 
tween two hard walls confining the polymers within a slit 
of width d. Using direct grand-canonical simulations of 
the full SAW polymer model, we computed the osmotic 
pressure exerted by the polymer-coils on the walls; the in- 
teraction free energy per unit area A, (3AF/A, is then ob- 
tained by integrating the osmotic pressure calculated for 
different values of the spacing d between the walls. These 
simulations are extremely computer intensive, and were 
only be carried out for L=100 |[L8| . In the "soft colloid" 
picture, the interactions of the polymer CM's with each 
other, @v(r), and with a wall, (3<fi(z), are calculated once 
with the HNC inversion procedures from the g(r) and 
p(z) of a full SAW polymer simulation at the bulk density. 
These are then used in grand-canonical MC simulations 
of soft particles between two walls, and in Fig. ^ they are 
compared to the 'exact" grand-canonical MC simulations 
of L=100 SAW polymers (for p/p* = 0.95) . The results 
are in good agreement, but the "soft colloid" calculations 
are at least two orders of magnitude faster. Contrary to 
the more widely studied case of colloid-colloid mixtures 
, the "exact" interaction exhibits no significant repul- 
sive barrier, whilst the "soft colloid" model leads to a flat 
maximum; the corresponding barrier height is, however, 
very small compared to the attractive minimum at con- 
tact, which agrees well with the "exact" data, as does 
the slope of the attraction. In fact, the repulsive barrier 
does not increase significantly with density |16[, and its 



origin can be traced to our use of the "potential overlap 
approximation" , namely that the interaction of the soft 
particles with two parallel walls a distance d apart can 
be written as the sum of the two individual wall-particle 
interactions. This is exact for simple liquids with true in- 
termolecular interactions, but not for polymers described 
by effective potentials, even if the polymers are ideal |L6| ]. 
For the sake of consistency, the MC simulations for the 
"soft colloid" model were carried out with effective wall- 
polymer and polymer-polymer potentials appropriate for 
L=100. However, we checked that the data obtained 
with effective interactions appropriate for longer poly- 
mers (L=500), which cannot be easily handled within 
the full segment model, are very close to the L=100 re- 
sults, so we are confident that we are close to the scaling 
regime. 
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FIG. 3. Depletion free-energy f3AF(d)/R 2 g for two 
plates separated by d. Circles are the "exact" MC simu- 
lations of SAW polymers, the diamonds denote MC sim- 
ulations of the "soft particles", the short dashed line 
denotes the wall-HNC approximation of Eq. (Q). The 
long-dashed and dash-dotted lines denote the AO approx- 
imations mentioned in the text. 

In Fig. [| we also compare two results derived in 
the spirit of the Asakura-Oosawa (AO) approxima- 
tion ^0|. The free-energy difference (3AF{z) is mod- 
eled by the density times the exact depiction volume, 
AVid(z), excluding one ideal Gaussian polymer of size 
R g , or by a popular phenomenological improvement [jp: 
f3AF(z)—f3Hi,AVid(z), where lit, is the bulk osmotic pres- 
sure of the interacting polymers. Note that for the den- 
sity under consideration, these approximations are seen 
to be very poor, both as regards the depth and the range 
of the depletion attraction. In fact, the range of the de- 
pletion interaction for interacting polymer coils is signifi- 
cantly reduced compared to the AO predictions, valid for 
ideal polymers. For low densities we find, as expected, 
that all the above approaches converge [|l6| . 

These observations can be understood within the "soft 
colloid" representation and the HNC approximation [Ml , 
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where the interacting free energy per unit area is given 
by: 

f3AF(z) 



= —p \ h(s)h(z — s)ds 



+ p h{z- s)[(34>{s)~ p(t> MF {s)]ds. (4) 



A 



Here h(z)=p(z)/p — 1 is the single wall density profile, 
I3(j){z) is the corresponding effective wall-polymer poten- 
tial, and f3(f> MF (z) is the corresponding potential of mean 
force. The first term on the r.h.s. of Eq. (0) is the 
density overlap approximation and would be the only- 
contribution in the case of ideal (Gaussian) polymer coils. 
The second term arises from the correlations between 
coils; this dominates the first term in the semi-dilute 
regime (p/p* > 1). The standard AO approximation 
|p0| may be derived from Eq. (||) by replacing the den- 
sity profile by a step-function of width R g in the first 
term of Eq. (Q) and neglecting the correlation term. In 
Fig. H we compare the HNC approach of Eq. (||) for the 
wall-wall interaction to the "exact" results and the MC 
simulations of the "soft colloids" . As was found for the 
homogeneous case and for the single wall, HNC works 
very well here, demonstrating that knowledge of (3v(r) 
and (3(f)(z) quickly leads to accurate predictions for the 
slit geometry, paving the way for the use of integral equa- 
tion techniques in other, more complex, geometries. 

To summarize, the coarse-grained representation of 
polymer coils as "soft" colloids has been shown to be very 
reliable, yielding pair distribution functions and concen- 
tration profiles which agree closely with the results for the 
full SAW segment model, while being much more efficient 
from a computational point of view. Much of the success 
of the coarse-graining lies in our finding that the "best" 
effective pair potential between CM's of neighboring coils 
does not depend strongly on polymer concentration, and 
is reasonably close to its p — > limit. Similar conclu- 
sions were reached in recent work on the phase-behavior 
of star-polymers, where the p — > limit of the pair po- 
tential was used to calculate the phase-behavior at finite 
concentration (22). Our results for the linear polymer 
case suggest that the full pair-potential for star polymers 
may not be strongly concentration dependent, and that 
our approach could be used for star polymers in confined 
geometries. 

Finally, we note that the "soft colloid" description 
is expected to work best in complex geometries where 
the curvature is not too large on the scale of R g , such 
as colloid-polymer mixtures where the colloid radius 
R < R g . For such systems, the "soft colloid" model 
may now be used in large scale simulations or fluid inte- 
gral equations of polymers in complex geometries, such 
as the structure Q, phase behavior M, interactions 
||, and metastability || of colloid polymer mixtures, 
which cannot be achieved with the detailed model of non- 
intersecting polymer chains. 
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